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Using a three-dimensional quasiclassical technique we explore molecular double ionization by a 
linearly polarized, infrared (800 nm) and ultrashort (6 fs) laser pulse. We first focus on intensities 
corresponding to the tunneling regime and identify the main ionization mechanisms in this regime. 
We devise a selection of observables, such as, the correlated momenta and the sum of the momenta 
parallel to the laser field as a function of the inter-electronic angle of escape where all the main 
mechanisms have distinct traces. Secondly, we address intensities above but close to the over-the- 
barrier intensity regime. We find a surprising anti-correlation of electron momenta similar to the 
experimental observations reported in Phys. Rev. Lett. 101, 053001 (2008). There, however, the 
anti-correlation was observed in very low intensities corresponding to the multiphoton regime. We 
discuss the mechanism responsible for the anti-parallel two-electron escape. 

PACS numbers: 33.80.Rv, 34.80.Gs, 42.50.Hz 



I. INTRODUCTION 

In the past two decades electron correlation has been 
established as the underlying mechanism for many im- 
portant phenomena arising from the interaction of strong 
laser pulses with matter. One of these phenomena is the 
dramatically enhanced multiple ionization yield of atoms 
(e.g. [lj]) and molecules (e.g. 0j) m intense laser pulses. 
Also called non-sequential multiple ionization, it is the 
laser driven recollision [!, H| of a field ionized electron 
with its parent ion that governs the multiple ionization 
process up to a certain threshold intensity Beyond 
the threshold intensity the contribution of recollision to 
multiple ionization vanishes and is replaced by a series 
of independent, sequential field ionizations. Using co- 
incidence imaging techniques such as COLTRIMS many 
experiments have succeeded in obtaining highly differen- 
tial kinematical details of electron correlation in the non- 
sequential intensity regime (e.g. [H, 0, H, H, H3, EH EH)- 
However, for higher intensities the available experimen- 
tal data is much less abundant. Correspondingly, the 
majority of theoretical work has concentrated on the non- 
sequential intensity regime (e.g. [H, El, EH)- 

Here, we report a classical study of electron correlation 
in N 2 molecules subjected to strong laser pulses at in- 
tensities well within the non-sequential double ionization 
regime (as indicated in Fig. [1]), as well as for intensities 
corresponding to the transition from pure tunneling to 
over-the-barrier ionization. 

First, we consider the intensity regime below the clas- 
sical over-the-barrier intensity. That is, one electron tun- 
nels in the field-lowered Coulomb potential, then accel- 
erates in the laser field and finally returns to the core to 
transfer energy to the remaining electron. In this inten- 
sity regime we find that the time of minimum electron 



approach is close to a zero of the field in agreement with 
the re-collision model Jj|. We differentiate the main en- 
ergy transfer mechanisms in the re-collision of the free 
electron with the parent ion using the time delay be- 
tween re-collision and ionization of the second electron. 
Furthermore, we devise a set of asymptotic observables 
where these mechanisms have distinct traces, such as the 
sum of the momenta parallel to the laser field as a func- 
tion of the inter-electronic angle of escape. 

Second, we consider intensities above, but close to, the 
over-the-barrier intensity threshold. This is an intermedi- 
ate regime where a transition takes place from correlated, 
or non-sequential, to independent, or sequential double 
ionization. We find that the instant of re-collision shifts 
to a maximum of the field in contrast to the common 
picture of non-sequential double ionization (NSDI) where 
the time of re-collision is close to a zero of the field. Fur- 
thermore, we find that the two electrons asymptotically 
depart in opposite directions. This is quite unexpected 
since for very high laser intensities one expects that both 
electrons ionize almost immediately parallel to the laser 
field and thus parallel to each other. We discuss how 
the escape of both electrons in opposite directions can 
be partly attributed to the laser field significantly in- 
teracting with both electrons before the re-collision - 
in contrast to smaller intensities where the effect of the 
field on the remaining electron can be neglected; another 
important factor is that the remaining electron is still 
bound and thus the effect of the Coulomb potential can 
not be neglected. 

Addressing the double ionization of strongly driven 
systems with first-principle quantum mechanical calcula- 
tions in all three spatial dimensions (3-d) is an immense 
task. Currently, 3-d ab-initio quantum mechanical calcu- 
lations are only available for the driven He atom [ijj . To 
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sions making it very important to accurately incorporate 
the Coulomb potential. 
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FIG. 1: Ratio of N2 double and single ionization yield as a 
function of laser intensity. Data points were experimentally 
determined in [l6| . At lower intensities non-sequential double 
ionization (NSDI) is the governing DI mechanism, though at 
higher intensities sequential double ionization (SDI) eventu- 
ally dominates. The threshold for over-the-barrier ionization 
(OTBI) of N 2 is about 2 x 10 14 W/cm 2 . 



cope with the highly complex task of tackling the double 
ionization of diatomic molecules many studies use numer- 
ical quantum approaches of reduced dimensionality (e.g. 
[H[). Others use judiciously chosen quantum models of 
reduced dimensionality (e.g. while some use semi- 

analytical quantum approaches in the framework of the 
so-called Strong-Field Approximation (SFA) (e.g. [2fj|). 
not fully accounting for the Coulomb singularity. 

We use a 3-d quasiclassical technique that we devel- 
oped for conservative systems, i.e. the sing le photon mul- 
tiple ionization of atomic systems pll. |22|. l23j|. Recently, 
we extended this technique to non-conservative systems 
to treat the correlated electron dynamics of the driven He 
atom 24j. Here, we further build on this work by tack- 
ling more than one atomic centers. Our method is nu- 



merically very efficient and treats the Coulomb singular- 
ity with no approximation in contrast to techniques that 
use "soft-core" potentials. Our method is also important 
for describing accurately effects such as the striking so- 
called "fin ger- like" structure which was recently observed 
in He [H 111 HI III- This structure was attributed to 
the strong interaction of the rescattering electron with 
the core — backscattering. In addition, accounting for 
the Coulomb singularity will be very important in pump- 
probe set-ups where VUV or XUV pulses are used. For 
these high frequency pulses, the excursion parameter of 
the electronic motion is smaller than the atomic dimcn- 



II. METHOD 

Our 3-d quasiclassical model entails the following 
steps: We first set-up the initial phase space distribu- 
tion of the two "active" electrons in the N2 diatomic 
molecule. Here, we consider only parallel alignment be- 
tween the molecular axis and the laser electric field. At 
intensities in the tunneling regime we assume that one 
electron tunnels through the field-lowered Coulomb po- 
tential. For the tunneling rate one can use quantum me- 
chanical or semiclassical formulas for diatomic molecules 
(see, e.g^[H, H3, HI HI)- We use the rate provided 
in ref. [29(. The longitudinal momentum is zero while 
the transverse one is provided by a Gaussian distribu- 
tion [3(j. This description is valid as long as the poten- 
tial barrier is not completely suppressed by the instan- 
taneous laser field E(t) = Eo(t) cos(wi). We consider 
the usual laser wavelength of 800 nm, corresponding to 
uj = 0.057a.u. (a.u. - atomic units). In our simula- 
tion the pulse envelope Eo(t) is defined as Eo(t) = Eq 
for < t < T and E (t) = E cos 2 {uj(t ~ T)/8) for 
T <t < 3T with T the period of the field. The thresh- 
old for over-the-barrier-ionization in neutral N2 with an 
ionization energy of I p \ = 0.5728 a.u. is reached at a 
field strength of E = 0.075a.u. (corresponding to roughly 
2 x 10 14 W/cm 2 ). 

Above 2 x 10 14 W/cm 2 the laser field allows an un- 
hindered electron escape and therefore the initial phase 
space is modeled by a double electron microcanonical dis- 
tribution [311 ]. However, it is important to note that in 
setting-up the initial phase space distribution we transi- 
tion smoothly from the tunneling to the over-the-barrier 
intensity regime. Namely, we assign a random number 
to the phase <j> of the laser field when the first electron 
is ionized, see [IH [32| . If the phase <j) corresponds to an 
instantaneous strength of the laser field E(4>) that leaves 
the electron below the barrier then we use the initial con- 
ditions dictated by the tunneling model. If the instan- 
taneous field strength pushes the barrier below the I p \ 
of that electron then we use the microcanonical distribu- 
tion to set-up the initial phase space distribution. This 
choice of initial conditions has proven successful in past 
studies (2o| in modeling the experimental ratio of double 
versus single ionization for long laser pulses [l(|. With 
our approach we ensure a smooth transition of the initial 
phase space distribution as we change the intensity. Even 
at an intensity of 3xl0 14 W/cm 2 still about 75% of the 
double ionization probability corresponds to trajectories 
initialized using the tunneling model, while 25% of the 
probability corresponds to trajectories initialized using 
the microcanonical distribution. 

After setting-up the initial phase space distribution we 
transform to a new system of coordinates, the so called 
"regularized" coordinates (33[. This transformation is 
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exact and explicitly eliminates the Coulomb singular- 
ity. This step is more challenging for molecular systems 
since one has to "regularize" with respect to more than 
one atomic centers versus one atomic center for atoms. 
We regularize using the global regularization scheme de- 
scribed in ref. .34]. Finally, we use the Classical Trajec- 
tory Monte Carlo (CTMC) method for the time propa- 
gation [35| . The propagation involves the 3-d four-body 
Hamiltonian in the laser field with "frozen" nuclei: 

2 2 1 1 

H = Vf ? il - 1 1 1 

^2 \R/2-n\ \-R/2- ri y 

+ 1 1 + (ri + r a ) • E(t), 

In - r 2 \ 

where E(t) is the laser electric field polarized along the 
z direction and further defined as detailed above, and R 
is the internuclear distance. 



III. RESULTS 

Re-collision as a subcycle, true attosecond process per- 
mits simulations using few cycle pulses to describe the 
main signatures of electron correlation that are mani- 
fested in multi-cycle pulse experiments Details of 
electron correlation as the branching ratio of the different 
NSDI pathways do depend on the pulse duration [3(| . In 
few cycle pulses the dominant channel of NSDI is the si- 
multaneous ejection (SE) of both electrons upon impact 
of the recolliding electron. The simultaneous ejection 
channel is responsible for the NSDI hallmark: a double 
hump structure in the sum of the momenta parallel to the 
laser field around Ay/Up. Here, 2^JTT V is the maximum 
velocity an electron can acquire from its interaction with 
the field and with U p = Eq/(Auo 2 ) being the ponderomo- 
tive energy. 

As the laser intensity is increased from 1 x 10 14 W/cm 2 
to 3 x 10 14 W/cm 2 a transition takes place as evidenced 
by the correlated momenta of the two escaping elec- 
trons shown in Fig. |2])a)-c). At 10 14 W/cm 2 both elec- 
trons escape with very similar momenta (Fig. l[2])a). Al- 
ready at 1.5 x 10 14 W/cm the distribution changes sig- 
nificantly and develops an X-like shape (Fig. (J2])b). How- 
ever, the inter-electronic angle of escape 9 (the angle 
between the asymptotic electron momenta) is generally 
less than 90°. In contrast, at 3 x 10 14 W/cm 2 the inter- 
electronic angle of escape becomes larger than 90°. In 
fact, the two electrons tend to escape with opposite mo- 
menta (Fig. ©c). This transition is also evident when 
one considers the distribution of 6 (Fig. ©c) medium 
panel) with a maximum around 165°. We note that our 
model does not account for depletion of the initial state 
in the-over-the-barrier regime. For that reason we only 
present the correlated momenta for intermediate inten- 
sities up to 3 x 10 14 W/cm 2 where the tunneling process 
still dominates as 3:1, thereby rendering the effect of de- 
pletion small. Finally, we note parenthetically that at 
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FIG. 2: Intensity dependence of electron correlation 
in the double ionization of N2 in a 3-cycle, 800 nm 
pulse. Columns a)-c) correspond to increasing intensities 
10 14 W/cm 2 , 1.5xl0 14 W/cm 2 , 3x 10 14 W/cm 2 , respectively. 
Top row: correlated momenta parallel to the polarization 
axis; Middle row: distribution of the inter-electronic angle 
of escape 8; Bottom row: sum of the parallel momenta vs. 
6. 

3 x 10 14 W/cm 2 over-the-barrier ionization of the first 
electron becomes possible, whereas NSDI is still the dom- 
inating mechanism in the liberation of the second elec- 
tron (compare Fig. [I}. 



A. Tunneling regime 

We now proceed to identify the double ionization path- 
ways for intensities in the tunneling regime. Two main 
mechanisms for non-sequential double ionization [l^ | are 
well established: simultaneous ejection (SE) upon rescat- 
tering and re-collision-induced excitation with subse- 
quent field ionization (RESI). 

In the SE mechanism both electrons are simultane- 
ously ejected with small but equal momenta, with the 
result that the asymptotic momentum components along 
the polarization axis are mainly dictated by the vector 
potential of the laser field at the time of ejection. This 
agrees with Fig. @a) where the maximum of the inter- 
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electronic angle of escape is around 30°. In RESI the 
rescattering electron only excites the remaining electron 
which is ionized at a subsequent maximum of the laser 

field [13, [II [11, Sit 

In Fig.[3]and Fig.@]we can clearly identify SE and RESI 
but also two other mechanisms. Namely, in Fig. ©c) we 
identify a NSDI pathway which starts out with excitation 
upon re-collision, but in contrast to RESI the excited 
electron ionizes around a subsequent zero of the field. We 
label this mechanism NSE2. In Fig. ((3j)d) we identify 
yet another mechanism that involves the formation of 
a doubly excited complex at the re-collision with both 
electrons ionizing at a later time as discussed in [li [3(| 
— we refer to this mechanism as the double excitation 
(DE) pathway. Previously, it was shown how the double 
ionization mechanisms have distinct traces in the spatial 
electron density [l9j . In the current work, we find that 
those different mechanisms also leave distinct traces in 
the observable momentum space. 

When electrons ionize through the SE pathway, they 
escape with large momenta and with inter-electronic an- 
gles less than 90°; in ionization through the RESI path- 
way the sum of the parallel components of the momenta 
is smaller and the electrons typically escape with angles 
larger than those for the SE mechanism. When the elec- 
trons escape through the NSE2 or the DE pathway the 
sum of the momenta is zero with large inter-electronic 
angles of escape. However, the two pathways have dif- 
ferent correlated momenta as shown in Fig. [3] and Fig. 0] 
(medium panel) — larger electron momenta are charac- 
teristic for the NSE2 mechanism. The SE and NSE2 are 
the dominant mechanisms for a 3-cycle laser pulse ac- 
counting for roughly 44% and 39%, respectively, for a 
laser intensity of 10 14 W/cm 2 and for 54% and 28%, re- 
spectively, for a laser intensity of 1.5 x 10 14 W/cm 2 , while 
RESI is around 10% in both cases. The contribution of 
the RESI and DE mechanisms increases significantly with 
increasing laser pulse durations. 

We note, that the identification of the different mech- 
anisms is based on the fact that the recolliding electron 
returns to the core, i.e. on the existence of an instant 
of closest approach to the remaining electron. We fur- 
ther verify the distinction between the different mecha- 
nisms by analyzing the time of ionization after re-collision 
shown in the bottom panels of Fig.[3]and Fig. [4] To deter- 
mine the time of ionization we make use of the "compen- 
sated energy" |4l[ for each electron. The compensated 
energy is the electron's potential energy with respect to 
the two nuclei plus its kinetic energy resulting from the 
canonical momentum due to the presence of the field. 
Here we define ionization as the instant when the com- 
pensated energy becomes positive and remains positive 
for all subsequent times. While being an approximate 
way to define the ionization time for each electron it 
proves useful for the identification of the different mech- 
anisms. 
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FIG. 3: For 10 14 W/cm 2 we identify four different double ion- 
ization mechanisms depending on each electron's time of ion- 
ization: the different columns refer to a) the SE, b) the RESI 
, c) the NSE2 and d) the DE mechanisms. Top row: the 
sum of the momenta of the two electrons as a function of 
the inter-electronic angle of escape; Middle row: the corre- 
lated momenta of the two electrons; and Bottom row: the 
ionization time of each electron in units of laser cycles. 
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FIG. 4: Same as for Fig. El but for 1.5 x 10 14 W/cm 2 . 

B. Over-the-barrier regime 

Now we are turning towards the results obtained at an 
intermediate intensity regime which is high enough that 
over-the-barrier ionization of the first electron becomes 
possible but low enough that the second electron cannot 
ionize without the energy transfer in a re-collision with 
the first electron (see Fig. ©c). The most striking fea- 
ture of this intensity regime is the strong anti-correlation 



5 



between the electron momenta, which resembles experi- 
mental observations by Liu et al. [Hj]. There, however, 
the anti-correlation was observed at intensities well be- 
low the classical threshold for re-collision double ioniza- 
tion. Since at these intensities a single re-collision does 
not provide sufficient energy to release a second electron 
a repeated energy transfer from the continuum electron 
via multiple field-assisted re-collisions was held responsi- 
ble for this effect. 

However, in the intermediate intensity regime cur- 
rently under consideration the energy of the recolliding 
electron easily exceeds the binding energy of the second 
electron. At 3 x 10 14 W/cm the impact energy at re- 
collision can be as high as 54 eV while the ionization 
energy of N^" is approximately I P 2 = 30eV. While this 
permits ionization also through shorter trajectories, as 
we discuss below, experiments have shown [10( that when 
excess energy is available only the minimum energy nec- 
essary for ionization is transferred in the collision and the 
recolliding electron retains most of the excess energy. 

To identify whether these "soft" collisions are indeed 
responsible for the observed anti-correlation we exam- 
ine in Fig. [3] the temporal evolution of the kinetic and 
potential energy of both electrons for laser intensities 
10 14 W/cm 2 and 3 x 10 14 W/cm 2 . The amplitude in the 
potential energy changes reveals that the collision for 
3 x 10 14 W/cm (thin lines) takes place further out from 
the nucleus if compared to the collision at 10 14 W/cm 
(thick lines). This interpretation is further corroborated 
by the presence of the "finger" -like structure in the cor- 
related momenta in Fig. (2)a) for 10 14 W/cm 2 and its ab- 
sence in Fig. (2)c) for 3 x 10 14 W/cm 2 . The "finger" -like 
structure was attributed in former studies of the strongly 
driven He atom to strong backscattering from the nucleus 
[13, Hj|. Finally, in Fig. |5] the temporal evolution 
of the kinetic energy shows that there is a smaller rel- 
ative change in the kinetic energy of the two electrons 
for 3 x 10 14 W/cm 2 around T/2 when compared to the 
kinetic energy change for 10 14 W/cm around 3/4T. 

At 3 x 10 14 W/cm 2 soft collisions with large impact pa- 
rameter seem to be favored due to the prevalence of short 
trajectories in the total double ionization yield. In Fig. [5] 
we identify the short trajectories by analyzing the tem- 
poral structure of the re-collision process. We plot the 
phase of the field when the first electron is "launched" 
as a function of the instant of double ionization (accord- 
ing to our previous definition of ionization time using the 
compensated energy). In the tunneling regime the recol- 
liding electron is "launched" shortly after the maximum 
of the laser field (Fig. ©a-b), i.e. the phase <f> of the 
laser field is close to zero (we use a cosine field). The 
majority of double ionization trajectories then recollide 
at around 3/4 into the laser cycle (zero of the laser field) 
corresponding to long trajectories. In contrast, in the 
over-the-barrier intensity regime the recolliding electron 
is "launched" when the phase of the laser field is larger, 
see Fig. ©c), corresponding to short trajectories. We 
find that the time delay between "launching" of the first 
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FIG. 5: Time dependence of kinetic energy (left panel) and 
potential energy (right panel) for the recolliding electron 
(black lines) and the initially bound electron (grey lines). 
Shown is the average over all double ionization trajectories 
at 10 14 W/cm 2 (thick lines) and 3 x 10 14 W/cm 2 (thin lines). 

electron (recolliding electron) and the first re-collision, 
i.e. the maximum of the electron-electron interaction en- 
ergy, decreases from 2/3T at I x I0 14 W/cm 2 to 1/2T at 
3 x 10 14 W/cm 2 . Although short trajectories have less 
than the maximum impact energy of 3.17U P they can be 
energetic enough at high intensities to release the second 
electron. However, the electron energy transfer becomes 
temporally less localized when the intensity is increased 
to 3 x 10 14 W/cm 2 . 

From the above it is clear that a transfer of energy 
through the re-collision event is still necessary for the 
bound electron to ionize for intensities in the interme- 
diate regime. However, the laser field for the higher in- 
tensities influences the motion of the bound electron even 
before the re-collision process. In contrast for the smaller 
intensities the asymptotic momenta components along 
the polarization axis are roughly dictated by the vector 
potential at the time of re-collision (around 2/3T). The 
significant role for higher intensities of the laser field even 
before the re-collision process can be also inferred from 
Fig. [7] Namely, depending on the time delay between the 
time of "launching" of the first electron and the time of 
double ionization, the two electrons escape either parallel 
or opposite to each other. In particular, the prevailing 
inter-electronic angle larger than 90° is due to the fact 
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FIG. 6: Instant of double ionization toi as a function of the 
phase cj> of the laser field when the first electron is "launched" . 
The laser intensity is a) 10 14 W/cm 2 , b) 1.5 x 10 14 W/cm 2 and 
c) 3 x 10 14 W/cm 2 . 



that the second electron ionizes approximately T/2 after 
the first electron is "launched" . With the vector poten- 
tial changing sign within half a period, this results in 
electrons escaping in opposite directions. 



Q, which would obscure the observed anti-correlation. 
However, at 3 x 10 14 W/cm 2 sequential double ionization 
can be neglected (compare Fig.[TJ. 

Further, our simulation does not include spatial inten- 
sity averaging as is virtually unavoidable in the experi- 
ment. Contribution of lower peak intensities to the over- 
all DI signal would partially mask the anti-correlation. 
This effect is expected to become increasingly important 
for higher intensities. 

Finally, one may legitimately ask to what extend the 
above finding is a molecular effect. In a recent study 
of the double ionization of the strongly driven He atom 
[2~4| the correlated momenta did not exhibit a pattern of 
electrons ejected in opposite directions. However, since 
our explanation of the anti-correlation does not necessi- 
tate molecular structure we would expect in principle a 
similar effect in atoms. Future studies in other atomic 
and molecular species can explore further whether the 
anti-correlation effect is generally observed for interme- 
diate intensities where there is an overlap of the over- 
the-barricr intensity regime and the regime where non- 
sequential ionization dominates over the sequential ion- 
ization. 



IV. CONCLUSION 
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FIG. 7: For 3 x 10 14 W/cm we show the distribution of the 
inter-electronic angle when the double ionization time is from 
T/4 to T/2 (circles); from T/2 to 3/4T (squares); from 3/4T 
to T (stars); and from T to 5/4T (diamonds). 



In conclusion, we have studied the non-sequential dou- 
ble ionization up to intensities where our classical model 
permits an intuitive definition of re-collision. In the 
tunneling regime we identify two new double ionization 
mechanisms and devise asymptotic observables with dis- 
tinct traces of these ionization mechanisms. For intensi- 
ties corresponding to the over-the-barrier ionization but 
still within the non-sequential double ionization regime 
we surprisingly find that both electrons are ejected with 
opposite momenta. In this intermediate intensity regime 
it is the field interaction with the second, bound elec- 
tron that permits soft re-collisions to assist the second 
ionization step. The combination of soft re-collision and 
timing between the two ionization steps forms the under- 
lying mechanism for the observed anti-correlation of the 
electron momenta. An anti-correlation pattern was also 
observed in 12j for small intensities where the mecha- 
nism responsible was identified as multiple inelastic field 
assisted recollisions. These results motivate a revision 
of the common notion that non-sequential double ioniza- 
tion necessarily means emission of both electrons with 
identical momenta. 



We will now discuss the limitations of our model with 
respect to the high intensity results. Our model does not 
incorporate sequential double ionization, i.e. tunneling of 
the second electron independent of the "launching" of the 
first electron. The independent tunneling of both elec- 
trons is known to yield uncorrelated electron momenta 
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